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ABSTRACT 



We present in this study a three step algorithm for the decomposition of 
arbitrary, three-dimensional, planar polygons into convex polygons. Through a 
series of translations and rotations, an arbitrary polygon is mapped onto the x-y 
plane, then broken into a set of convex polygons, and finally mapped back to the 
polygon's original coordinate system for filling and display by special graphics 
hardware. 



3 



TABLE OF CONTENTS 



I. INTRODUCTION 6 

II. BACKGROUND 7 

A. DEFINITIONS 7 

B. DECOMPOSITION 7 

C. RESEARCH FACILITIES 9 

III. STEP I: MAPPING THE 3D POLYGON TO THE 2D X-Y PLANE 12 

A. THE NORMAL TO A POLYGON 12 

B. ANGLES ABOUT THE AXES 14 

C. TRANSFORMATION MATRIX 16 

IV. STEP II: POLYGON DECOMPOSITION 18 

A. CLOCKWISE ORDERING 18 

B. VERTEX EXAMINATION 19 

1. Test One 19 

2. Test Two 19 

3. Test Three 22 

C. SELECTION COMPLETION 24 

V. STEP III: MAPPING TO THE ORIGINAL COORDINATE SYSTEM .. 26 

VI. IRIS IMPLEMENTATION 27 

A. MATRIX IMPLEMENTATION 27 

B. DOUBLY LINKED LIST 27 

C. VERTICES STACK ! 28 

VII. LIMITATIONS AND CONCLUSIONS 30 

LIST OF REFERENCES 33 



4 



APPENDIX - ALGORITHM SOURCE CODE 34 

INITIAL DISTRIBUTION LIST 76 



5 



I. INTRODUCTION 



Computer hardware advancements applied to the field of computer graphics 
have recently made possible the real-time display and animation of three- 
dimensional object models. In most currently produced graphics systems, the 
three-dimensional (3D) object models are comprised of sets of 3D polygons. For 
years the key bottleneck to generating real-time images of 3D models has been the 
limited speed with which the 3D polygons comprising the model have been filled. 
This bottleneck has been alleviated recently with the advent of special hardware 
designed to fill polygons at rates exceeding 44 million pixels per second [Ref. l] . 
For real-time animation of 3D object models, such polygon fill rates are essential. 

The majority of the graphics systems that provide polygon fill hardware only 
provide hardware for filling convex polygons. Concave polygons passed to such 
hardware usually are filled incorrectly, if at all. To utilize the capabilities of such 
hardware for the general polygon case, that is, a mix of concave and convex 
polygons, it is necessary to have some means of decomposing an arbitrary 3D 
planar polygon into a set of 3D planar convex polygons completely covering the 
area defined by the original polygon. The following study attempts to provide 
such an algorithm. We then discuss the implementation and performance 
characteristics of that algorithm on the Silicon Graphics, Inc. IRIS workstation. 



6 



II. BACKGROUND 



A. DEFINITIONS 

A polygon is defined as a geometric figure that is described by an array of x. 
y. and z coordinates for its vertices. A two dimensional polygon lying in the x-v 
plane is a special case of a three dimensional polygon in which the z coordinates 
are all zero. A convex polygon is one that has no points of negative curvature 
along its boundary [Ref. 2: p. 217]. This means that a line can be drawn from 
any point inside the polygon to any of its vertices and have the line lie entirely 
within the polygon boundary [Ref. 3: p. 350]. A concave polygon fails the above 
test. i.e. we cannot draw a line from every point inside the polygon to every 
vertex and have the line remain within the polygon boundary. However, any 
point inside a concave polygon can be shown to lie inside a convex polygon 
created by some of the vertices and sides of the original polygon [Ref. 3: p. 330]. 
Figure 2.1 provides examples of convex and concave polygons. 

B. DECOMPOSITION 

To fill a concave polygon utilizing hardware that only correctly fills convex 
polygons, we must first break down the polygon into simpler convex polygons that 
together cover the area of the original. The process of decomposing a concave 
polygon does not necessarily yield a single unique set of convex polygons, but the 
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Figure 2.2 - Two Equivalent Decompositions 
of a Concave Polygon. 
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result of two different decompositions of the same polygon should be equivalent as 



shown in Figure 2.2. 

During the operation of decomposing a polygon, there are two main checks 
that must be conducted to ensure complete and proper decomposition takes place. 
The first check is to ensure that no area of a created convex polygon falls outside 
the boundary of the original concave polygon. The second check is to ensure that 
the concave polygon is completely covered by the total areas of the created convex 
polygons. Figure 2.3 provides examples of the two checks. 

The polygon decomposition algorithm produced for this study is comprised of 
three basic steps. The first maps the three dimensional coordinates of the polygon 
into a two dimensional plane. This process generates a polygon that is more 
easily examined and decomposed. The second step conducts the actual 
decomposition an performs the two decomposition checks. The third step maps 
the coordinates of the created convex polygons back to the the original orientation 
in three-space. At this point, the created convex polygons are then ready for 
passing to the graphics system’s polygon fill hardware. 

C. RESEARCH FACILITIES 

The IRIS (Integrated Raster Imaging System) 2400 Graphics Workstation, 
manufactured by Silicon Graphics, Inc., is the target system for our polygon 
decomposition algorithm. By incorporating custom VLSI chips to replace slower 
graphics software, the IRIS system can accept polygons in user-defined 



9 




original polygon, vertices could be chosen that 
yield an incorrect decomposition. 




Figure 2.3b - The second decomposition check ensures 
the concave polygon is completely covered by the 
decomposed polygons. In this case all vertices 
and edges were utilized, however a hole is still 



left. 



10 



coordinates, and transform them to screen coordinates for display on a color 



graphics monitor that provides 1024 x 768 resolution. This use of special 
hardware yields processing speeds that are some 100 times faster than similar 
operations carried out in software. 
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III. STEP I: MAPPING THE 3D POLYGON TO THE 2D X-Y PLANE 



A. THE NORMAL TO A POLYGON 

The first step of the polygon decomposition algorithm, is comprised of several 
small steps. The goal of these steps is to map the polygon to the x-y plane. First, 
three non-colinear points are selected from the array of x, y. and z coordinates of 
the polygon. These are used to define the equation of the plane for the 3D 
polygon, and the normal to that plane (Figure 3.1). 

Ax + By + Cz = D 
N - A_» + Bj + Cj 

The coefficients A. B, and C are calculated from the coordinates of the three 
selected points by the following equations: 



a = (y 0 - y t )(z 2 - zJ - (y* - yi)(zo - *i); 



B = (z 0 - Zi)(x 2 - x,) - (z 2 - zjfxo - Xj); 



C = (xo - x,)(y 2 - y t ) - (x 2 - x^yo - y,). 



The coefficients are then used to solve for D of the equation of the plane. 
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Three vertices 0, 1, and 4 

are selected to determine the 
Normal to the plane. 



Figure 3.1 The Normal to a Polygon 
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B. ANGLES ABOUT THE AXES 



Once the equation of the plane is determined, each point in the polygon input 
array is checked to ensure that the polygon is in fact planar. If the polygon 
passes this test, then the angles between the normal and the axes of the 
coordinate system are found. From these angles, the necessary rotations, about 
the x and y axes, to bring the polygon into the x-y plane, are calculated. The 
following equations [Ref. 4: pp. 24-59] are used to calculate the angles shown in 
Figure 3.2a. 



cos a 



A 

v/A' + B* + C a 



cos f) = 




C* 



cos 7 VFTFTc 5 

For the polygon to lie in the x-y plane, the angle between the normal and the 
z axis is either zero or 180 degrees. Necessarily then rotations are conducted to 
bring the angles between the normal and the x and y axes to 90 degrees each. 

First a projection of the normal onto the x-z plane is made as shown in Figure 
3.2b. The dot product of that projected normal N' and the z axis is calculated. 
The dot product is divided by the product of the magnitudes of the two vectors 
yielding the cosine of the angle of rotation about the y axis (<&). 
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N 




Figure 3.2 Angles about the Axes 
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cos ^ = NT' 



2 

Itf'l |Z| 



A second projection of the normal onto the y-z plane is made and the angle of 
rotation about the x axis (0) is calculated in the same manner as above (Figure 
3.2c). 



C. TRANSFORMATION MATRIX 

To actually rotate the polygon into the x-y plane, a transformation matrix 
must be created that is used in the multiplication of the polygon's coordinates. 
We first translate the polygon to the origin by adding the negative x, v, and z 
coordinates of one of the polygon's vertices. Assuming 1, m, and n are the 
coordinates selected, the matrix becomes: 



10 0 0 
0 10 0 
0 0 10 
—1 — m — n 1 

A rotation matrix for the angle <z> about the y axis is then created 



cos^ 0 sin^ 0 
0 10 0 
— sin^ 0 cos+ 0 
0 0 0 1 

and concatenated with the translation matrix. This is followed by a multiplication 
with a rotation matrix for the angle 6 about the x axis. 



16 



10 0 0 
0 costf — sinO 0 
0 sinfl cos 0 0 

0 0 0 1 

The result of these matrix concatenations is the following transformation matrix 
that is used to map the array of vertices onto the x-y plane. 



COS^ 


— sin^sinfl 


— sin^cosfl 


0 


0 


cos0 


sin0 


0 


sin^ 


— cos^sinfl 


cos^costf 


0 



— lcos^— ncos^ — mcos0+(lsin^— ncos^)(— sin0) — msintf+(lsin^— ncos^)(cosfl) 1 

Once we have this matrix set up. call it M loplane , we then use it to transform 
the polygon’s coordinates via a simple row vector by matrix multiplication. The 
result of this multiplication is a transformed set of coordinates, in which the z 
coordinates are zero. The x and v values of this transformed array are then used 
in the decomposition step of our algorithm. We note at this point that for the 
mapping of the coordinates back to the original vertices (Step III of the Alg.), 
that the index into the transformed array is also the index for the untransformed 
point in the original array of vertices. 
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IV. STEP II: POLYGON DECOMPOSITION 



The decomposition step divides the polygon into multiple convex polygons. 
In so doing, the algorithm performs the two checks to ensure that no area outside 
the polygon is shaded and that the created convex polygons completely cover the 
area of the concave polygon. This step requires that the vertices of the polygon 
be ordered in a clockwise fashion, with respect to the +z directed normal to the 
x-y plane. Thus a determination of the ordering must be made and if necessary 
the array of vertices reversed. 

A. CLOCKWISE ORDERING 

To determine if the vertices of the polygon are ordered clockwise, the 
polygon's area is calculated using the formula [Ref. 5: p. 369] below. (Note r,,y, 
are the polygon's 2D coordinates.) 

y((xoyi) + (xiy 2 )+ *“ +( x n-iyn)+(^nyo)-(yoXi)-(y 1 x 2 )- ••• -(yn-i*J-(yn*o)) 

The above equation yields a negative area for clockwise and a positive area for 
counterclockwise polygons. If the polygon is ordered counterclockwise then the 
array of x, y, and z coordinates that define the polygon is reversed before the 
selection procedure is begun. 
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B. VERTEX EXAMINATION 



The selection of vertices for a convex polygon is conducted by examining each 
vertex of the concave polygon with respect to the previously selected vertices of 
the developing convex polygon. This is begun by initially selecting two adjacent 
vertices, then in clockwise order, examining each of the polygon’s remaining 
vertices to build as large a convex polygon as will fit within the boundary of the 
concave polygon. 

A series of tests are conducted in the selection of vertices for the convex 
polygon. As vertices are selected, the line segments connecting these become the 
edges of the convex polygon and define its boundary. 

1. Test One 

The first test of a vertex is if it lies to the left or right of each line 
segment of the developing convex polygon. The vertices of the concave polygon 
are examined in clockwise order. Thus if a vertex lies to the left of any of the line 
segments of the developing polygon, then a concave angle is necessary to include 
the vertex in the polygon. As shown in Figure 4.1a, vertex 2 lies to the left of the 
line segment 0-1. A concave angle ( p in Figure 4.1b) is necessary to add vertex 2 
to the developing polygon. Therefore vertex 2 is not selected and the process 
moves on to perform the same test on vertex 3. 

2. Test Two 

A vertex that passes the first test may still require a concave angle to be 
included in the developing convex polygon. In Figure 4.2. the vertices 0, 1, 2, 3, 
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Figure 4.1 - (a) Vertex 2 lies to the left of the 
line segment 0-1, (b) Concave angle ( p ) would 

be required to connect Vertex 2 to the devel- 
oping polygon, (c) Vertex 2 is skipped and 
the examination is continued with vertex 3. 
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1 



2 




Figure 4.2 - (a) Vertex 5 lies to the right of each 
line segment of the developing polygon, (b) Concave 
angle ( p ) is required to connect Vertex 5 to Vertex 
0, (c) Concave angle (p) can be identified since 0 

lies to the left of line segment 4-5, (d) Vertex 4 

is skipped and the test repeated with line segment 
3-5, (e) Vertex 3 is skipped and the test repeated 

with line segment 2-5, (f) Result of the tests are 

that vertices 0, 1, 2 and 5 are selected. 
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and 4 have been previously selected for a convex polygon with vertex 5 under 
examination. Vertex 5 passes the first test since it lies to the right, or inside of 
each of the line segments. However, a concave angle (/> in Figure 4.2b) is again 
required to connect vertex 5 with vertex 0 of the developing polygon. A concave 
angle of this type is identified, as shown in Figure 4.2c, if the initial vertex of the 
developing convex polygon lies to the left of the line segment created by the 
vertex under examination and the last selected vertex of the developing polygon. 
In Figure 4.2c. the initial vertex. 0, lies to the left of line segment 4-5. To remedy 
this, the last selected vertex of the developing convex polygon, vertex 4, is 
removed and the test is again conducted with the line segment formed by the 
vertex under examination and the last selected vertex (Figure 4. 2d A: 4.2e). This 
process is repeated, removing vertices from the developing polygon, until the 
initial vertex no longer lies to the left of a new line segment. 

3. Test Three 

The third test conducted is to ensure that no line segment created in the 
convex polygon intersects any of the edges of the concave polygon. Such an 
occurrence results in area outside the concave polygon being included in the 
convex polygon and an inaccurate polygon fill taking place. In Figure 4.3, the 
previously selected vertices are 0, 1, 2, and 3. Vertex 7 is under examination and 
passes the two earlier tests. However, the selection of vertex 7 results in the 
convex polygon 0, 1, 2, 3, 7 which contains area outside the concave polygon. A 
situation such as this is identified since line segment 3-7 intersects line segment 
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Figure 4.3 - Tests must be conducted for the 
intersection of developed line segments with 
the edges of the original concave polygon. 
Otherwise Vertex 7 would be an acceptable 
vertex for inclusion in the developing 
polygon 0, 1, 2, 3... 
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13-0 of the concave polygon. When this occurs, the vertex under examination is 
removed from the convex polygon and the next vertex of the concave polygon is 
examined. 

C. SELECTION COMPLETION 

The selection process for the developing convex polygon concludes when all 
vertices of the concave polygon have been examined. The convex polygon is then 
ready for the third step of our algorithm, mapping to the original polygon’s 3D 
coordinates. Figure 4.4 is an example of a concave polygon that was decomposed 
to the three shaded convex polygons using the above three tests. We have a 
"hole" left in the concave polygon even though all vertices and edges of the 
concave polygon have been utilized. Holes such as this are prevented by taking 
the initial and last vertices of a developed convex polygon as the starting vertices 
of a convex polygon to be developed. In Figure 4.4, vertices 0 and 2 of polygon 0, 
1, 2 are taken to develop convex polygon 0. 2. 4 to fill in the "hole". Once this is 
accomplished, step two of the algorithm is repeated by selecting a vertex of the 
concave polygon that has not yet been utilized in a convex polygon. The 
decomposition process is complete when there are no vertices or edges of the 
concave polygon that have not been utilized in one or more convex polygons. 
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1 3 




Figure 4.4 - All vertices and edges of 
the polygon have been used, yet a hole 
remains . 
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V. STEP III: MAPPING TO THE ORIGINAL COORDINATE SYSTEM 



The second and third steps of the decomposition process are interleaved, i.e., 
each created convex polygon is immediately mapped back to the original 
coordinate system. This is accomplished by maintaining the original array of 
vertices that the user provided for the three-dimensional polygon. As a convex 
polygon is produced, it is actually represented by the indices of its vertices. Actual 
vertex coordinates are extracted from the vertex array, as required, using the 
index. In this manner, the transformed coordinates can be utilized during 
decomposition and the index can be used to obtain the original coordinates for 
final display. By this method, there is no consequential need to translate and 
rotate the vertices back to the original coordinate system. 
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VI. IRIS IMPLEMENTATION 



Implementation of the decomposition algorithm on the IRIS system is via 
approximately 1600 lines of C code (Appendix). Several data structures and 
various IRIS system hardware features are utilized as described below. 

A. MATRIX IMPLEMENTATION 

In the development of the transformation matrix, IRIS system calls for 
translation and rotation are used to alter the top of the system matrix stack. 
These routines perform the matrix concatenations needed to translate and rotate 
the polygon to the x-y plane. To convert the polygon vertices to 2D coordinates, 
four vertices at a time are loaded into a matrix. The x, y, and z coordinates go 
into the first three columns with l’s filling the fourth column to provide the 
necessary homogeneous coordinates. This is placed on the system’s matrix stack 
and multiplied with the transformation matrix using the IRIS’s Geometry 
Engines. The transformed coordinates are removed from the matrix stack and the 
process is repeated until all vertices have been transformed. 

B. DOUBLY LINKED LIST 

A data structure for each vertex of the polygon is maintained in our system. 
This structure contains the vertex’s index into the vertex array. As actual 
coordinates are needed, the vertex’s index is used to extract the x, y, and z 
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coordinates from this array. A "used” flag indicates the vertex's use in a 



developed convex polygon, with a slope and y-intercept stored for the line 
segment the vertex forms with its predecessor in the array. A forward pointer 
connects the structure with the next vertex in clockwise order, while a back 
pointer connects the structure with its preceding vertex. This produces a doubly 
linked circular list of vertex structures that permits forward clockwise traversal of 
the polygon and rearward counterclockwise traversal (Figure 6.1). 

C. VERTICES STACK 

To manage the examination and selection of vertices for a developing convex 
polygon, a stack is utilized. The size of the vertices stack is determined at run 
time when contiguous memory is dynamically allocated, based on the number of 
vertices in the concave polygon. As a vertex is examined, its index is placed on 
top of the stack along with the slope and y-intercept of the line segment it creates 
with the vertex below it. If the vertex passes all selection tests, it remains on the 
stack and a new vertex is pushed on top. On the other hand, when a vertex fails 
a selection test, it is popped off the stack and the next vertex to be examined is 
placed on top. Once the selection of vertices for a convex polygon is completed, 
the vertices are passed to the IRIS system hardware and the stack is cleared to 
continue with the decomposition process. 
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(a) Doubly linked circular list of structures, 



Rearward 

Pointer 


Index 


Used 

Flag 


Slope 



Y-intercept 



Forward 

Pointer 



(b) Vertex structure 



Figure 6 . 1 
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VII. LIMITATIONS AND CONCLUSIONS 



This study has presented an algorithm for the decomposition of an arbitrary, 
planar polygon into a set of convex polygons, thus eliminating the need to 
manually perform such decomposition. With the algorithm’s dynamic allocation 
of memory, the size polygon that can be filled is limited only by the amount of 
memory available for use. The algorithm indirectly provides an extended 
capability over the IRIS system's limit of 384 vertices in one polygon (Ref. 6: p. 
2 ). 

However due to the time required to test each vertex for concavity, the 
algorithm is not suggested for use where real-time response is important. To 
demonstrate the overhead involved in testing for concavity, a convex polygon was 
filled and displayed 500 times by the algorithm. The same polygon was then 
filled utilizing only the IRIS's polygon fill hardware. The results are provided in 
Figure 7.1 for various polygon sizes. 

Due to the three primary steps necessary for decomposition and the tests 
conducted on each vertex, the speed of the algorithm depends entirely upon the 
number of vertices and concave angles of the polygon. As such, the success of the 
algorithm's use in real-time display and animation is based upon a balance 
between the complexity of the concave polygons used and the speed requirements. 
We suggest that the created concave polygons be computed once and saved. 



30 



Seconds 




Figure 7.1 - Time required to fill and 
display 500 convex polygons with 
the number of vertices indicated 
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For IRIS usage, this is perhaps best accomplished through the use of the "object" 
mechanism of the IRIS graphics library. 
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APPENDIX - ALGORITHM SOURCE CODE 



**************** 4 - ********************************************* 

globals.h 

****** 4 -****************************************x*************** 

/ 



# include < gl.h > 

f include <device.h> 

# include <stdio.h> 

# include <math.h> 

/************************************************************** 

GLOBAL CONSTANTS 

/ 

£ define NULL 0 
f define BOTTOM 0 

# define new(x) ix *) malloc(sizeof(x)) 

GLOBAL TYPES 

/ 



typedef int boolean; 

;’ + Define a structure to contain data about a vertex of the polygon. 

These structures will be used to create a linked list. * 

typedef struct point { 
struct point *back; 
int index; 
float slope; 

Coord intercept; 
int direction; 
boolean failed; 
boolean segment used; 
struct point ^forward; 

}; 

typedef struct point point; 
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/* Define a structure for storing data about each vertex. This will 
be used to create a stack for later use. * 

typedef struct stack { 
int index; 
float slope; 

Coord intercept: 
int direction; 

}; 



typedef struct stack Stack; 

/***% ******** X******** **************** ************************* 

GLOBAL DATA STRUCTURES AND VARIABLES 

*************************************************************** 

Stack *stack; 

point *start ptr, ^current ptr; 
int TOP; 

boolean clockwise_error_flag; 
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I * ************** ******************************* ****** ********** 

ROUTINE: convert. c 

*X**X*****************XX**********±x**x***x:jk**x.**.*.*-******.*.***** f 

This routine takes th initial array of vertices, passed by 
the user, and makes the calls to other routines to perform the 
necessary translation and rotations to transform the array into 
the x-v-0 plane This new array is then used by the decomoser to 
breakdown the polygon into convex polygons. 



cpolf(total pts. pts) 

int total pts: 

Coord ptsj] 3 ; 

{ 

int ptl, pt2, pt3. cnt; 

float A. B. C, D, anglex, angley. anglez. *neu pts: 

Matrix transform 1; 

Matrix points: 

stack=(Stack *) calloc(total_pts -h 1. sizeof(Stack)): 

t %***********%******** *%**?£** ******* ***%*%**^%* *%*%%*** 

dynamically allocate memory to create an array to hold the 
vertices of the newly transformed polygon. 

** ************************* -XL*****-#.*** **************** ********** j 



new pts= (float *) calloc (3*total pts, sizeof(float)); 

I ************************************************************** 

get 3 non-colinear vertices from the polygon to be used for 
calculating the normal to the polygon. 



extract 3 pts(total pts, pts, &ptl, &pt2, &pt3): 

I ************************************************************** 

calculate the coefficients for the equations to the plane 

of the polygon and the normal to that plane. 

*************************************************************** 

determine coeffs(total pts, pts, ptl, pt2, pt3, &A, &B, &C, &D); 
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check to ensure all vertices of the polygon lie in the same 
plane. If so then calculate the necessary angles of rotation 



if(planar polygon(total_pts, pts, A. B. C, D)) { 

calculate the rotation about the v axis necessary to 
bring the normal into the v-z plane. 



angles from normal(A, B, C, 5 y\ &anglex, Bangley); 

create the transformation matrix necessary to multiply 
with the array of vertices to generate the rotation. 



build transform matrix(transforml, angley, ’y\ pts, ptl); 

perform the matrix multiplication with the array of 
vertices and store the new coordinates in the new array 
new pts. 



map pts(total pts, pts, transform!, points, new pts); 



now extract three non-colinear vertices from the new 
array of vertices, new pts. 



extract_3_pts(total pts, new pts, &ptl, &pt2, &pt3); 

/XXXXX*XXXXXXXXXXXXXXXXX*XXXXXXXXXXXXXXX*XXXXXXXXXXXXXXXXXXXXXX 

since the polygon has been rotated about the y axis we 
calculate the coefficients for the equations to the 
polygon one more time. 

determine_coeffs(total_pts, new pts, ptl, pt2, pt3, &A, &B, &C, &D); 
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calculate the rotation about the x axis 
angles from normal(A. B. C, ’x\ fcanglex, &angley); 



k^********** 



create the transformation matrix 

* * + i 



build transform matrix(transforml. anglex, ’x\ new_pts. ptl); 

/ * * * > 

perform the matrix multiplication 

* ***a 



map_pts(total pts, new pts. transform!, points, new pts); 



' * +c *********:» 



k*** ******** :M 



call the routine which will perform the decomposition 
upon the newly transformed polygon which now lies in the 
x-y-0 plane. (thus to be handled as a 2D polygon) 



k******?M 



k * * * * -i 



split(total pts, new' pts, pts); 



once the dynamically allocated memory for the array 
new-pts has served its purpose, the memory is freed. 



cfree(new r pts. 3*total pts. sizeof(float)); 
cfree(stack, total pts, sizeof(Stack)); 



} 



if not planar then the decomposition is not attempted. 



return; 
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/************************************************************** 

ROUTINE: extract. c 



This routine searches the array of vertices to find three points which 
do not lie on the same line and returns the indices of these points 



*******: 4 :******************************************************* 



/ 



extract 3 pts(total pts, pts, ptU pt2. pt3) 

int total pts; 

Coord pts[] [3] ; 
int *ptl, *pt2, *pt3; 

{ 



int count, finished: 

*ptl= 0; *pt2= 0; *pt3 = 0: counts 0: finished^ 0; 

/* The first vertex in the array (index=0) is selected for the first point. */ 
while(count < total pts &&; [finished) { 

/* We then search the rest of the array of vertices to locate two more */ 
count-b= 1; 
if(*pt2 == 0) { 

/* If the second point has not been selected yet then we check that 
x, y. and z coordinates, of the vertex currently under exam, are 
different. */ 

if(pts[*ptl O' != ptsscount] fOj || 
pts: *ptl] [l] != ptsjcountjjl] || 
pts[*ptl][2] != ptsicount][2] ) 

/* If so then the vertex becomes the second point. */ 

*pt2= count; } 

else if(*pt3 == 0) { 

/* If the third point has not yet been selected, then we check that 
the coordinates of the vertex under examination are different 
from both the first and second point. */ 

if((pts[*ptl][0 != ptsicount][0] 
pts *ptl l] != pts'count] 1 || 
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pts[*ptl 2 != ptsicount] 2 ) 

(pts *pt2, 0 != ptsicount 0 
pts[*pt2] 1 != pts count 1 
pts[*pt2][2 != pts count] 2 )) 

/* If the coordinates are different then the vertex is selected 
as the third point. */ 

*pt3= count; } 



else 

finished= 1; 



} 



return; 



} 
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/*******************************************************¥****** 

ROUTINE: is planar.c 

*************************************************************** / 

This routine takes each vertex of the polygon and uses it to solve the 
equation for the plane as calculated earlier. If every vertex lies 

in the plane, then a value of 1 is returned. 

*************************************************************** / 

planar polygon( total _pts, pts. A, B. C, D) 

int total pts; 

Coord pts[] [3] ; 
float A. B. C, D; 

{ 



int count, planar: 

count— 0; planar— 1: 

while(count < total pts &:& planar) { 

/* The equation for the plane is setup to equal zero if the coordinates 
of the vertex lies in the plane. A range of .02 on either size of 
zero is allowed to handle round off error. This test is done for 
each vertex in the polygon. If any vertex fails the test then the 
function returns (0) for non-planar. */ 

if((A*(pts count 0 )+ B*(ptsicount T )+ C*(pts count 2j)- D) > .02 j 
(A*(pts|count 0])+ B*(pts count][l )-b C*(pts count] [2])- D) < -.02) 
planar= 0: 

count-h= 1: 



} 

return(planar); 



} 
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/I************************************************************** 



ROUTINE: coeffs.c 

This routine uses the three vertices selected from the polygon to 
calculate the coefficients that will be used in both the equation 
of the plane the polygon lies in and the normal to that plane. 



determine coeffs(total pts, pts, ptl, pt2, pt3, A. B, C, D) 

int total pts: 

Coord pts 3'; 
int ptl, pt2, pt3: 
float * A. *B, *C, *D; 

{ 

float xO. xl, x2, yO, yl, v2. zO. zl, z2:; 

x0= pts ptl 0 : 
y0= pts ptl [1 ; 
z0= pts[pt l] [2] ; 

xl= pts pt2 [0 ; 

>1= ptsjpt2][l;; 
zl = ptsi pt2 [2j ; 

x2= pts pt3 0 : 
y2= pts pt3 1 ; 
z2= pts pt3] [2] ; 

*A= (yO-v 1 )*(z2-zl) - (y2-y l)*(z0-zl); 

*B= (z0-zl)*(x2-xl) - (z2-zl)*(x0-xl); 

*C= (x0-xl)*(y2-y 1) - (x2-xl)*(y0-y 1); 

*D= (*A)*pts ptl jOj + ( * B) *pts[pt 1] jl] + (*C)*pts[ptl‘ 2j; 

return; 



} 
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/************************************************************** 



ROUTINE: angles.c 

This routine utilizes the coefficients, determined from the normal to 
the polygon, to determine the angles of rotation necessary about the x 
and y axes. These rotations are necessary to bring the polygon into the 
x-y plane. This routine is called the first time to calculate the 
rotation about the y axis. The second time it is called is to determine 
the x axis rotation. The parameter "axis" is used to pass which angle 
is to be calculated. 



angles from_normal( A, B. C, axis, anglex, angley) 

float A. B. C: 
char axis; 

float *anglex, * angley: 

{ 



float pi, degx. degy, degz; 
pi= 3.1415926536; 

degx= acos(A/sqrt( A* A^B*B+C*C))*360/ (2*pi): 
degy= acos(B/sqrt(A*A+B*B^C*C))*360/ (2*pi); 
degz= acos(C/sqrt(A^A-B*B-C*C))*360/(2 + pi); 

if(A==0 && B==0) { 

/* simply translate the polygon to the origin */ 

*anglex= 0.0; 

*angley= 0.0; } 

else if(A==0 && C==0) { 

/* rotate polygon 90 degrees about x axis */ 

*anglex= 90.0; 

*angley= 0.0; } 

else if(B==0 && C==0) { 

/* rotate polygon 90 degrees about y axis */ 

*anglex= 0.0; 

*angley = 90.0; } 

else { 

/* must calculate amount of rotation about the x & y axes */ 
4 anglex= acos(C/sqrt(B*B+C*C))*360/(2*pi); 

* angley = acos(C/sqrt(A*A^C*C))*360/ (2*pi); 



43 



if(axis == V) { 
if(degx < 90.0) 

*anglev = -(*angley); 
*anglex= 0.0; } 
else { 

if(degy > 90.0) { 
*anglex= -(*anglex); 
*anglev= 0.0; } 

} 



/* convert the angles into tenths of degrees as needed by the IRIS system. */ 
*anglex= 10*( + anglex); 

*angley= 10*(*angley); 

return; 



} 
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ROUTINE: transmatx.c 

This routine takes a matrix structure and initializes it to an identity 
matrix. Then the necessary translation values and rotation values are 
added to the matrix to create the resultant transformation matrix. 

ft *?************************************************************ i 



build transform matrix(transforml, angle, axis, pts, ptl) 

Matrix transform 1: 
float angle; 
char axis; 

Coord pts 1 3 ; 
int ptl; 

j 

\ 

int i, j; 

/* create the identity matrix * 
for(i=0: i<4; i-r = l) { 
for(j=0: j<4: j+=l) { 
if(i==j) 

transforml i] [j = 1; 
else 

transforml^ 1 j = 0; 

} 

} 

pushmatrixQ; 

loadmatrix(transforml); 

/* add in the rotation values */ 
rotate((int)angle, axis); 

/* add in the translation values */ 
translate(-pts[ptl 0j, -pts(ptl][l], -pts[ptl 2]); 

get matrix (transforml); 
popmatrix(); 

return; 

} 
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ROUTINE: mapover.c 

*************************:«*:***************** 4 :****************** / 

This routine places up to 4 vertices into a matrix at one time then 
performs a matrix multiplication with the transformation matrix to 
produce x, y . and z coordinates for the polygon that is translated 
and rotated some amount about the x and y axes. These new coordinates 

are restored in the array new pts for further manipulation. 

*************************************************************** 



map pt total pts, pts, transforml, points, new' pts) 

int total pts: 

Coord pts 3 ; 

Matrix transforml, points: 
float new_pts' 3 ; 



{ 



int count. count2, cnt; 
count= 0: count2= 0: 
while(count < total_pts) { 
cnt= 0; 

while(cnt < 4) { 

/* place the coordinates of a vertex into the matrix */ 

pointsjcntl 0 = pts count 0 ; 

poin tsjcnt] [ 1] = pts count lj; 

pointslcnt] |2 = pts count 2]; 

pointscnt] ;3 = 1.0: 

if((count4-= 1) < total pts) 
cnt— = 1; 
else 
break; 



} 



/* save the state of the matrix stack and perform the multiplication */ 
pushmatrixQ; 



loadmatrix(transforml); 
multmatrix (points); 
get matrix (points): 
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popmatrix(); 
cnt= 0; 

while(cnt < 4) { 

/* move the newly calculated coordinates back into the vertices array */ 
new pts c.ount2] 0 = points cnt][0j; 
new ptsicount2] lj— points cnt] i l]; 
new pts count2],2 = points cnt] [2^ 

if((count2— = 1) < count) 
cnt-f= 1; 
else 
break; 



} 



} 



return; 



} 
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ROUTINE: split, c 

This routine tak^s the array for the transformed polygon and calls 
the necessary fu lions and routines to decompose the polygon 
and pass the vertices of the convex polygons to the IRIS system 
for filling. 

*******¥***:*:*:fc************************************************* / 



split(total_pts. pts. old_pts) 

int total pts: 

Coord pts 3j; 

Coord old _pts[] 3 ; 

{ 



Object decompose polygonQ: 

Object whole_thing: 
boolean clockwise: 
point *pl, *p2: 

/*^****^****************^***********^************************** 

Create the doubly linked circular list, "main chain”, for the 
concave polygon represented by the pts array. Also initialize 
flags and ptrs that will be used in decomposing the polygon. 

******4:*^****4:^*********************5tc************************** / 

/ 

clockwise= clockwise ordering(total pts, pts); 

create main chain(clockwise, total pts, pts); 

/ *************************************************** *********** 

Decompose the main polygon and build the set of simple convex 
polygons. 

*************************************************************** f 



decompose_polygon(total_pts, pts, old pts); 
return; 



} 
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/************************************************************** 

ROUTINE: main_chain.c 

***************************************************************. 

Create the doubly linked circular list, ’’main chain’ 1 , for the 
concave polygon represented bv the pts array. Also initialize 
flags and ptrs that will be used in decomposing the polygon. 

create main chain(clockwise. total pts. pts) 

boolean clockwise: 
int total pts; 

Coord pts [^ 3;; 



{ 



int countl: 
point *pl, *p2: 

/* allocate memory for a structure to contain data on the first 
vertes of the polygon * 

pl= new(point); 

/* set pointers to this structure and initialize all values */ 

start ptr= pi; 
current ptr= pi; 
pl->index= 0; 
pl->direction= 0; 
pl->failed= FALSE; 
pl->segment _used= FALSE: 

/* if polygon vertices are ordered clockwise then we index the vertices 
array beginning at 0 and incrementing from there. Other we must begin 
at the other end of the vertices array to extract our coordinates. In 
this way a forward move along our doubly linked list will equal a 
clockwise move around the polygon. */ 

if(clockwise) 
countl= 1; 
else 

countl= total pts- 1; 

/* for each vertex of the polygon we build a point structure, calculate 
the slope, intercept, initilize flags and connect the structure 
via pointers to its preceeding and succeeding vertices in the 
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polygon. The result is a doubly linked list of structures. * 
do{ 

p‘2= new(point); 

Pi ->forward= p2; 
p2->back= pi: 

Pl= P 2: 

pl->index= countl: 

determine slope(pts pl->back->index 0 , ptS(pl-> back->index 1, 
ptsipl->index 0 . pts pl->indexj 1 , &(pl->slope), 

<Mpi- > intercept). &(p Indirection)); 
plnfailed= FALSE; 
plnsegment used= FALSE; 
if(clockwise) 
countl^= 1: 
else 

count 1 -= 1: 

}while(countl< total pts && countl > 0); 

/* connect the first vertex structure to the last vertex of the polygon 
and calculate the slope and intercept of the line segment the two 
create. * 

Pi nforward= start ptr; 
start ptr-> back= pi: 
pl= start ptr; 

determine slope(pts pl«> back- > index 0 , ptsjplnbacknindex] [l , 
pts'plnindex 0 . pts plnindexjjlj, &(plnslope). 

&( P i -> intercept). &(p Indirection)); 



return; 



} 
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/ ******************************************************* jt****** 

ROUTINE: clock wis.c 

This routine calculates the area of the arbitrary planar 
polygon once it has been transformed to the 2D x-y plane. 

A negative area indicates a clockwise ordering of the vertices 
while a positive area indicatesa a negative orderig of vertices 
A "0” is returned for counterclockwise ordering while a ”1” is 
returned for clockwise ordering. 



clockwise ordering(total pts, pts) 

int total pts; 

Coord pts[] 3]; 

{ 



float area: 
int i. result; 

area= 0.0: 

for(i^0; i< total pts- 1; i-h= 1) { 
area= area + pts i 0 * pts [i-hl] [ 1 ; 

} 



area= area -r pts total pts- 1 0 * pts 0 1 : 

for(i=0; i< total pts- 1; i-h= 1) { 
area= area - pts [i i 1 * pts i iH- 1 [0 ; 

} 



area= area - pts; total pts- 1 1 1 * pts[0] [0] ; 

if(area >= 0) 
result= 0; 
else 

result= 1; 
return (result); 



} 
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/************************************************************** 



ROUTINE: decompose, c 

*************************************************************** / 

Decompose the main polygon and as convex polygons are determined 
these are passed through calls to the IRIS system to be filled. 

if:***:************:**********:*****************:*******:*****:******:** / 



Object decompose polygon(total pts. pts. old pts) 

int total pts: 

Coord pts[] 3j; 

Coord old ptsj !3 : 



{ 



boolean no point found, search incomplete, adjust needed. OK; 
boolean skipflag. skip_search: 
int ptl, pt2, poly count, value; 
point ^working ptr: 

clear_stack(); 
get_pt(total_pts, pts); 
get pt(total pts, pts); 

search incomplete= TRUE; 
skip search — TRUE: 
poly count= 0: 



/**************a;*********************************************** 

Once an initial two points are taken from the concave polygon, 
the process moves through the main chain to select out points 
for constructing the simple convex polygons. Calls are made to 
various routines to check for validity of pts that are accepted 
into a polygon. 

*************************************************************** / 



while(search_incomplete) { 
adjust needed — TRUE; 

OK- TRUE; 
while(OK) { 
get pt(total pts, pts); 

/* we first check to see if we have made a complete loop of the 
vertices. */ 

if(loop completedQ) { 
working ptr= start ptr; 
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/* find the structure in the linked list that matches the 
vertex on top of the stack * 

while(working ptr->index != stack TOP 1 . index) { 
working ptr= working_j>tr->forward; 

} 



/* check to see if vertex creates a line segment with the 
vertex below it. on the stack, that intersects any edges 
of the concave polygon * , 

if(stack'TOP-l .index != working ptr->back->index) { 
value= intersecting lines( pts i stack ' TOP I. index; Oj , 

pts stack TOP .indexl 1 , pts stack TOP-1 .index] Oj, 
ptslstackTOP-1. index 1 , stack TOP .slope, 
stack TOP .intercept, total pts, pts); 

} 

if(value) { 

/* if vertex creates an intersecting line segment then we 
remove it from the stack and reset the pointer to the 
linked list * 

current _ptr= start ptr; 

while(current ptr->index != stack[TOP-l .index) { 
current ptr= current ptr->forward: 

} 



current ptr= current ptr->forward; 
releaseptQ; 
release pt(); } 

else { 

/* since the line segment create is OK, we simply remove the 
vertex from the stack since it is the same as the one on 
the bottom of the stack, namely the beginning vertex */ 

OK= FALSE; 
release_pt(); 

} 

} 

/* if we have not completed the loop then we must conduct the 
tests to ensure valid selection of vertices for a convex 
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polygon */ 

else if(pt_out of_bounds(totai_pts, pts)) { 
release pt(); 

} 

} 

/* if at least three points have been selected then... */ 
if(valid pts()) { 

/* save the beginning and ending vertices for future reference */ 

pt 1 = stack BOTTOM!. index; 
pt2= stack TOP .index: 

/* pass the vertices to IRIS system hardware for filling and display * 

define_polygon(old pts); 
poly count — = 1: 

* if this is a normal decomposition and not a 11 fill in” polygon to 
cover holes left in the concave polygon the the do the following * 

if(skip_search) { 

/* set the flags that indicate next time through this loop we 
will be working on a "fill in" polygon */ 

adjust needed= FALSE; 
skip search= FALSE; 

/* adjust the pointers to the linked list so next time through 
the loop we have the start and end vertices of the previous 
polygon to begin working with */ 

working ptr= start_ptr; 
while(working_ptr->index != pt2) { 
working_ptr= working _ptr-> forward; 

} 



/* if the start and end vertices of the previous polygon are 
adjacent on the concave polygon then forget about building 
a "fill in" polygon */ 

if(working ptr->forward->index != ptl) { 
while(current ptr->index != ptl) { 
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current _ptr= current _ptr->forward: 

} 

get pt(total pts, pts); 
current ptr= working ptr; 
get pt(total pts, pts); } 
else { 

adjust needed= TRUE; 

} ' 

} 

} 

/* if the polygon is not fully decomposed yet we select 

two new vertices to continue the decomposition process */ 

if(adjust needed) { 
skip_search= TRUE; 
if(object complete( total pts)) 
search incomplete— FALSE; 
else { 

adjust search area(pts); 
clear stackQ; 
get pt( total pts, pts); 
get_pt(total pts, pts); 

} 

} 

} 



return; 



} 



/**********************************:k**' ****:**:*:f************4:*** 

ROUTINE: clearstack.c 

***********************:«********:*****************************:** / 

/ 

This routine simply resets the stack pointer to -1. In this 
manner when the stack has its first item pushed on, the pointer 
will accurately reflect a value of 0. (BOTTOM OF STACK) 

********************************************************** 4 :**** / 



clear stack() 

{ 

TOP= -1: 

} 
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/************************************************************** 

ROUTINE: get pt.c 

*************************************************************** / 

This retrieves the next point in the main chain, as identified 
by the current ptr. and places the point on the top of the stack 
In doing so it increments the TOP variable. Also as the point 
is retrieved, the slope, intercept, and direction of the line 
segment it creates with the vertex below it on the stack is 
calculated. 

*********************************************** **************** / 



get pt (total pts, pts) 

int total_pts; 

Coord pts[]|3 ; 

{ 

TOP— = 1: 

if(TOP <= total pts) { 

* select the vertex in the linked list that is currently pointed to * 
stack TOP .index= current ptr->index; 

/* and then advance the list pointer * 
current ptr = current ptr- > forward; 
if(TOP~> 0) { 

/* take the x. v , and z coordinates for both the vertex on top of 
the stack and the vertex just below it. to calculate the slope 
and v-intercept of the line segment they create */ 
determine_slope(pts stack TOP-1 .index, [0 . ptsjstacktTOP-1. index] [l , 
pts stack TOP .index][0 ,, ptsjstackjTOP]. index] [l], 
fcstackTOP .slope. &stack TOP .intercept, 

&:stack TOP .direction); 

} 

} 

else { 

printf( M STACK overflow occured while retrieving another point"); 



return; 



} 
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/it:*************************************:**-******* ^i************** 

ROUTINE: releasept.c 

X***%********************X****X*****************X*****X******** 

This routine removes the item on top of the stack by 
decrementing the top of stack pointer TOP. 



release pt() 

{ 



if(TOP >= BOTTOM) 

TOP-= 1: 

else { 

printf( ?, STACK underflow occured while trying to release a non- M ); 
printf("existence point. "): 

} 

return: 

} 
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ROUTINE: slopes. c 

±*X***X********************************-********x*************** 

This works through the pts array, and the stack array. By 
taking value of the TOP of stack and the point beneath it. 
the corresponding coordinates are obtained from the pts array 
and a slope, intercept and direction are calculated. These 
values are then stored in the stack. 

****************:*:*********** ******************:*: +*** 4 -* ********** 



determine slope(xl, vl. x2, y2. slope, intercept, direction) 

Coord xl, vl. x2, v2: 
float *slope: 

Coord ^intercept: 
int ^direction: 



{ 



Coord dx, dy: 

/* calculate the delta x and delta y for the to points * 



dx= x2- xl; 

dy= y2- yl; 



/* if the line segment is vertical then use the value 99999.0 to 
indicate such (avoiding division by zero problems) * 

if(dx == 0) { 

*slope= 99999.0; 

*intercept= x2; } 
else { 

*slope= dy/dx: 

*intercept= v2 - (dy/dx) * x2: 

} 



/* based upon the changes in delta x and delta y we determine if we are 
traversing down a slope or up a slope */ 

if(dy < 0) { 
if(dx < 0) 

*direction= 1; 
else 

^direction— 0; 

} 
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else { 

if(dx > 0) 

* direction= 0; 
else 

+ direction= 1: 

} 



return; 



} 
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/I****************************************************** ******:** 



ROUTINE: compare. c 

****************************************************+*:********* 

This routine determines if a point lies to the left or right of 
a line segment. This is done bv using the x and y coordinates 
of the point with the slope, intercept, and direction of the 
line segment to create a parallel line segment. The v intercept 
(or in the case of vertical line, the x intercept) is used in 
conjunction with the direction of traversal along the line 
segment to find which side of the line segment, the point lies 
on. 



compare boundry (xl.y 1. slope, intercept, direction) 

Coord xl. yl: 
float slope: 

Coord intercept: 
int direction: 



{ 



float test_value; 
int value; 

/* if the line segment is vertical then we use the x-intercept to 
determine if a vertex lies to the left or right of another line 
segment */ 

if(slope == 99999.0) 
test_value= xl; 
else 

test_value= yl- slope * xl; 

/* based upon traversal up or down a sloping line we determine if 
a vertex is to the left or right of the line segment */ 

switch (direction) { 

case 0: if(test value > intercept) 
value= 1: 
else 

valuer 0; 
break; 

case 1: if(test_value < intercept) 
value= 1: 
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else 

value= 0; 
break; 

} 

return (value): 



} 
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/************************************************************** 



ROUTINE: ptOOB.c 

*************************************************************** 

This routine takes the vertex index from the top of the stack 
and uses its x. y coordinates to determine if the point lies 
out of bounds with respect to the convex polygon being 
created. 

*********************************************** X'**** *********** 



pt out of bounds(total_pts, pts) 

int total pts; 

Coord pts 3 ; 

{ 

Coord xl, vl: 
int ptr. value; 
boolean not finished: 
point * working ptr: 

* obtain the coordinates for the vertex on top of the stack * 

xl = pts stack TOP .index 0 ; 
vl = pts stacklTOP . index 1; 
ptr = TOP- 1: 
value= 0; 

* use the above coordinates to determine if the vertex lies to 
the left (out of bounds) with respect to any line segment 
created thus far in the selection of vertices */ 

while (ptr > BOTTOM lvalue) { 
value= compare_boundry(xl,y l,stack[ptr .slope. stack ptr . intercept. 

stack [ptr 1 . direction); 
ptr-— 1: 

} 



/************************************************************** 

If point is not out of bounds with any line segments created 
then we check to see if the segment produced by this new point 
addition causes the start point for the polygon to fall out of 
bounds. 

*************************************************************** / 



if(lvalue) { 

xl= pts[stack[BOTTOM].indexJ 0 ; 
yl= ptsjstackiBOTTOMj. index 1 ; 
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value= compare boundry (xl,y 1, stack TOP .slope.stackiTOP]. intercept, 
stack [TOP .direction): 
if(value) { 
release pt(); 

current ptr= current ptr->back: 
working ptr= start ptr; 

} 

} 



If the newly created line segment does not cause the initial pt 
of the developing polygon to fall out of bounds, then we check 
to see if the line segment intersects any line segments of the 
concave polygon. 



if(!value) { 

working ptr= start ptr; 

while(working ptr->index != stack TOP .index) { 
working ptr= working ptr->forward; 

} 

if(stack TOP-1 .index != working ptr->back->index) { 

* we use the coordinates for the vertex at the top of stack 
and the vertex below it to create a line segment which we 
then check for intersection with any of the edges of the 
original concave polygon * 

value= intersecting lines(pts stack TOP .index] 0 . 

pts stack TOPi.index][l , ptslstack[TOP-l .index [0 . 
pts stack TOP-1 .index, [1], stack TOP .slope, 
stack TOP .intercept, total pts, pts); 

} 

} 



If the newly created line segment does not intersect any line 
segments of the concave polygon then we check to see if the 
line segment of the original polygon, that the top of stack 
point heads, causes the point below top of stack to fall out. 

if(!value) { 

if(current_ptr->back->back->index !— stack[TOP-lj. index) { 
xl = ptsjstack TOP-1 !. index] 0]: 
yl= pts stack TOP-l]. index 1 ; 

valuer compare boundry (xl,y 1, current ptr- > back- > slope. 
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} 



} 



current ptr->back->intercept, 
current ptr->back-> direction); 



If the point is determined to be out of bounds, then the 
appropriate pointers are adjusted to continue the decomposition 
process. Otherwise the pt is a valid addition to the stack 

'/ 



if(value) { 

working ptr= start ptr; 
not _finished= TRUE: 
do { 

if(working_ptr->index == stack'TOP .index) { 
not finished= FALSE: 
if(working_ptr->failed == FALSE) 
working ptr->failed = TRUE;} 
else { 

working ptr= working ptr- > forward: 

} 

} while(not finished): 

} 

return(value); 



} 
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ROUTINE: loop complet.c 

X****X***X*********X************±**>±*x**X *********************** 

This routine checks to see if the vertices of the concave 
polygon have been completely traversed while developing a single 
convex polygon. This is done when the vertex on top of the 
stack is the same vertex on the bottom of the stack. When this 
occurs a value of 1 is returned. 



int loop completedQ 

{ 

int value; 
if(TOP> -1) { 

if(stack TOP .index == stack BOTTOM .index) { 
value= 1; } 
else { 
value= 0: 

} 

} 

else { 

printf(”Error in stack index while checking TOP & BOTTOM values”) 
value= -1; 

} 

return (value): 

} 
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/**$*%** ************************************************* ****** 



ROUTINE: adjust. c 

This routine searches the circular list of vertices, checking 
the various condition flags, to locate a vertex or line segment 
of the originial polygon which has not been utilized. As one is 
located, the list’s currency pointer is moved to that structure 
in the list so that further polygon decomposition may continue 



adj ust search _area( pts) 
Coord pts[] 3 : 

{ 



boolean done: 
point *pl: 

done= 0: 
pl= startptr: 

do { 

if(pl-> failed == 1) 

/* check to see if the line segment to either side of the vertex 
has been used in a previously developed polygon. x / 

*f(!(pl -> forward-> segment used) !(pl-> segment used)) { 
current ptr= pi: 
pl-> failed= 0; 
done= 1; } 
else { 

pl= pl-> forward: 

} 

else { 

pl= pl-> forward: 

} 

} while(!done && pi != start ptr); 
return: 



} 
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.'************************************************************** 

ROUTINE: valid _pts.c 

*************************************************************** 

This routine checks to ensure that at least 3 vertices have 
been determined for the convex polygon. If not then the 
value of 0 is returned. 

************* ************************************************** / 



boolean valid pts() 

{ 



int value: 



if(TOP >= 2) 
value= 1: 
else 

value= 0: 
return (value); 



} 
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/******************************¥*** * 4 - ** + + 4 -***** + **************} 

ROUTINE: definepolv.c 

**************************************************************4 

This routine takes the vertices accumulated on the stack and 
uses the IRIS features pmv and pdr to actually perform the 
display of the polygon created. The index for each vertex 
is used to access the original array of points to get the x. y 
and z coordinates used. 

:i.4::4;K*4c:*:4:*:********4 : ****4- - *4:*4-'***4- + ******4:**:f¥:»:*4:**4-4-4 : -**:’*:*4:4:4:4:4:*** 



define_ polygon (pts) 
Coord pts 3 : 

{ 



int count: 
point *track ptr: 

track ptr= start ptr; 

/* pass the starting vertex for the newly created convex 

polygon to the IRIS hardware with a point move command * 

pmv(ptsistack BOTTOM .index) 0 . 
pts stack BOTTOM .index' 1 . 
pts stack BOTTOM .index 2]); 

/* locate this vertex in the linked list * , 

while(track_ptr->index != stack [BOTTOM .index) { 
track_ptr= track_ptr->forward: 

} 



/* set the flag to indicate the vertex has been used 
in a convex polygon */ 

if(track _ptr->back->index == stack TOP .index) 
track_ptr->segment_used = TRUE: 

/* for each vertex of the convex polygon, pass its 
coordinates to IRIS hardware with the point draw 
command */ 

counts 1; 

while(count <= TOP) { 
while(track ptr->index != stackicount). index) { 
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track ptr= track ptr->forward: 

} 

if(track ptr->back->index == stack count- 1 .index) 
track_ptr->segment_used= TRUE: 
pdr(pts stack count .index] 0 . 
pts stack [count .index] [l . 
pts stack count .index. 2]): 
count^= 1: 

} 

/* pass the command to the IRIS hardware that the 
fill the completed convex polygon * 

pclosQ; 

clear stackQ; 

return; 



} 
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ROUTINE: obj complet.c 

This routine checks each vertex and edge of the original 
concave polygon to see if they have all been utilized in the 
process of decomposition. If so a value of 1 is returned. 

**%.***# t ****** %* 4-** **********%* ********x*%* *%*$■-************** ** 



int object complete(total pts) 
int total pts; 



{ 



point *pl: 

int count, value; 

count= 0: 
value= 1: 
pl= start ptr: 

/* search through the linked list of structures, for the 
vertices, and return the index of the first vertex that 
creates an edge of the concave polygon not yet utilized * 
do { 

if( P i -> segment used) { 
pl = pl-> forward; 
count -t-= 1; } 
else { 

value= 0: 

} 

} while(count < total pts &.& value == 1); 
return (value); 



} 
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/****************************************************-*******:**»: 

ROUTINE: intersecting. c 

*;x*:*************************************************sc*****:**:^:*^ 

This routine takes two vertices and the slope and intercept 
of the line segment they create to determine if it intersects 
any of the edges of the original polygon. If so then a value 
of 1 is returned. 

*^f****** ***********************************************: *;■**:*:*•**;*•. 



intersecting lines(x2, y2. xl, yl. slope* intercept* total pts, pts) 

Coord x2. v2, xl* yl: 
float slope: 

Coord intercept: 
int total pts: 

Coord pts ( 3; 



{ 



point * pi: 

float x intersect, y intersect. x_min. x_max. y_min* y max, fudge: 
int value: 

* define a value to offset round off error * 

fudge= 0.000001; 

value= 0; 
pi = start ptr; 

for each edge of the concave polygon we check if either vertex 
lies on the edge. If not then we determine the intersection 
of the line* defined bv the edge, with the line defined by the 
two vertices. Finally we determine if this intersection lies on 
both line segments. If so then we have found an actual intersection. */ 

do { 

/* check if first vertex lies on the polygon's edge being examined */ 

if( ( (xl ! = pts|pl->indexj[0]) || (yl != pts[p 1- > index] [ 1 ) ) && 

((xl != pts|pl->back->index O’) || (yl != pts[pl->back->index li))) { 

/* check if second vertex lies on the edge */ 

if(((x2 != ptsipl->index 0 ) |l (v2 != pts j p 1- > index] [ 1 ] ) ) && 

((x2 != pts|pl->back->index 0 ) (y2 != pts, pi- >back->index 1 ))) { 
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/* if the lines are not parallel then ... 
if(slope 1= pl->slope) { 

/4 ' if the line created by the vertices is vertical 
if(slope == 99999.0) { 

/* determine the coordinates of the intersection 
x intersect^ x2: 

y intersect= pl->slope * x intersect-!- pl->intercept; } 

/* if the edge of the polygon is vertical * 

else if(pl->slope == 99999.0) { 

/* determine the coordinates of the intersection * 
x intersect^ pts pl->index 0 : 
y intersect= slope * x intersect— intercept; } 

/* if neither line is vertical 4 

else { 

/* determine the coordinates of the intersection 4 
x_intersect= (intercept - (pl->intercept))/((pl->slope) - slope); 
v intersect= slope * x intersect^- intercept; 

}” 

/* determine the x range that the intersection must occur in 
to lie on the line segment created by the two vertices *7 

max_min(x2.xl.&x_max,&x_min); 

/* check if the intersection falls within this range */ 

if((x intersect >= x_min-fudge)&&(x_intersect <= xrnax+fudge)) { 

/* since the intersection is within this x range we now determine 
the x range for the edge of the concave polygon */ 

max_min(pts[pl->index] 0], pts pl->back->index f0., &x max. &x_min) 

/* check if the intersection falls within this range */ 

if((x_intersect >= x_min-fudge)&&(x_intersect <= x_max+fudge)) { 

/* now that the intersection falls within the x range of both 
line segments, we must calculate the y range for each 
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segment and perform the checks 



/* y range for the vertices line segment * 
max min(y2,yl.&y max,&y min); 

if((y intersect >= y min-fudge)<§c&(y intersect <= v max^fudge)) { 

/* y range for the concave polygon edge * 

max min(pts pl->index 1 ,pts pl->back->index 1 ,&y max.fc y min) 
if((y intersect >= y min-fudge)&&(y intersect <= y max— fudge)) { 
value= 1; 

/* the intersection exists for the actual line segments and 
therefore the vertex on top of the stack must be removed */ 

} 

} 

} 

} 

} 

} 

} 

pl= pl->forward; 

} while((pl != start_ptr)&&:(value == 0)); 
return (value): 



} 
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ROUTINE 



max min.c 




Given two float values, pi and p2. this routine determines 
which is larger of the two. 



max_min(pl, p2. max, min) 
Coord pi, p2, "max. *min: 



if(pl < p2) { 

*max= p2: 
*min= pi; } 
else { 

*max= pi: 

* min= p2: 

} 

return; 




{ 



} 
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